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The conductance in an extended multiband Hubbard model describing linear arrays of up to 
ten quantum dots is calculated via a Lanczos technique. A pronounced suppression of certain reso- 
nant conductance peaks in an applied magnetic field due to a density-dependent spin-polarization 
transition is predicted to be a clear signature of a coherent "molecular" wavefunction in the array. 
A many-body enhancement of localization is predicted to give rise to a giant magnetoconductance 
effect in systems with magnetic scattering. 

PACS numbers: 73.20.Dx, 71.27. +a, 73.40. Gk, 78.66.-w 
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Arrays of coupled quantum dots Jy^ can be thought of 
as systems of artificial atoms separated by tunable tun- 
nel barriers. Two complementary theoretical approaches 
have been useful in describing such systems in the limit 
H where charging effects are important but interdot 
tunneling is incoherent and in the limit H of coher- 
ent ballistic transport, with charging effects neglected. 
However, recent improvements in fabrication and exper- 
imental techniques should make it possible to probe a 
third regime, where both interaction and coherence ef- 
fects play nontrivial roles. In this regime, the system of 
coupled quantum dots behaves like an artificial molecule, 
and must be described by a coherent many-body wave- 
function HH]. In this Letter, we describe some strik- 
ing characteristic signatures of such a coherent molecular 
wavefunction in the low-temperature magnetotransport 
through an array of quantum dots. Our theoretical pre- 
dictions should be experimentally testable in currently 
available GaAs quantum dot systems. 

An important consequence of coherent interdot tun- 
neling is the formation of interdot spin-spin correlations 
|jj analogous to those in a chemical bond at an energy 
scale ~ t 2 /U, where U = e 2 /C g is the charging energy of 
a quantum dot and t = (h 2 /2m*) J d 3 x **(x)V 2 *j(x) 
is the interdot hopping matrix element, <IVj being elec- 
tronic orbitals on nearest-neighbor dots. In an applied 
magnetic field, the Zeeman splitting drastically modi- 
fies the many-body wavefunction of the array when the 
Zeeman energy g^sB ~ t 2 /U, and we find that certain 
resonant conductance peaks are suppressed (or, in some 
cases, enhanced) by several orders of magnitude com- 
pared to their size at B = 0. For quantum dots elec- 
trostatically defined in a 2D electron gas , no signifi- 
cant modification of the wavefunction of a single quantum 
dot would occur for magnetic fields of this magnitude in 
the plane of the dots (the case we consider), so that the 
standard Coulomb blockade-based transport theory ||] 
which neglects coherent interdot tunneling would predict 
no interesting magnetic field dependence of the conduc- 
tance. Such spin-spin correlations are intrinsic many- 
body effects which are non-perturbative in the Coulomb 



interaction, and can not be explained by ballistic trans- 
port theories 0] either. We therefore believe that ob- 
servation of the predicted dramatic magnetic field effect 
on low-temperature transport through coupled quantum 
dots would represent the clearest possible signature of 
the formation of an "artificial molecule." 

The system we wish to model consists of a linear ar- 
ray of quantum dots electrostatically defined |[]|^] in a 
2D electron gas, with a magnetic field in the plane of 
the dots. We neglect intradot correlations (a reason- 
able approximation if the number of electrons per dot is 
not too small HH) and focus instead on collective phe- 
nomena in the array. The electron-electron interactions 
in the array are described by a capacitance matrix 
Cij : we assume constant capacitances C g between each 
quantum dot and the macroscopic metallic gate which 
defines its confinement potential, and a capacitance Cj, 
which is a function of gate voltage and may include im- 
portant quantum mechanical corrections |9|,|ifj||, between 
nearest neighbor quantum dots. The electronic orbitals 
in the confining potential of an isolated dot are taken to 
be nondegenerate with level spacing A, and the hopping 
matrix element t„ between (nearly) degenerate orbitals 
on nearest-neighbor dots is assumed largest, all others 
being neglected. The Hamiltonian is 



H- 



na-Cjno- + H.C. 



where Cj n<7 is the creation operator for an electron of spin 

a in the nth orbital of the jth dot, hj = J2 n a ^jno^jna, 
and the sums run from n = to M — 1 (the M orbitals 
nearest the Fermi energy Ep), i,j = 1 to L, and a = ±1. 
In the strongly-correlated regime, the interaction term in 
Eq. (Q) cannot be treated perturbatively. We therefore 
employ a Lanczos technique |ll] to compute the exact 
many body ground states of (Jl) for arrays of 5 to 12 
quantum dots with 1 to 5 electronic orbitals per dot. 
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The array is coupled to noninteracting leads via a tun- 
neling Hamiltonian with matrix elements t£ oc t n which 
couple electrons in the nth orbital of the 1st (Lth) dot to 
the right (left) lead. The capacitance to the leads is ne- 
glected. In the limit AE > k B T > HT, where AE is the 
energy level spacing in the array and T is the tunneling 
rate of electrons out of the array, the linear response con- 
ductance is determined by ground state to ground state 
transitions, and is given by 
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{ON\t r /c\ {L)n jO N ^)\ 2 p r ' l (E° N - 

i)/H, p r ' l (e) being the density of states in the 
leads. For the case B > 0, the ground state is non- 
degenerate and Am{p) = —f'{E^ — E^_ 1 ), while for 
B = 0, the ground state is spin-degenerate when N is 
odd and A N (p) = 2/[k B T('S + 2e~ XN + e XN % where 
x N = -E% + El^/keT. Eq. (§) is derived 

by the method of Refs. ||§. 

Fig. ^ shows the conductance through a linear array 
of 10 quantum dots with Ci = as a function of the 
chemical potential \x in the leads, whose value relative 
to the energy of the array is controlled by the gate volt- 
ages. The two Coulomb blockade peaks in Fig. [ljare split 
into multiplcts of 10 by interdot tunneling, as discussed 
in Refs. ^,0]. We refer to these multiplets as Hubbard 
minibands. The energy gap between multiplcts is caused 
by collective Coulomb blockade ||, and is analogous to 
the energy gap in a Mott insulator @]. The heights of 
the resonant conductance peaks in Fig. |l](a) can be un- 
derstood as follows: Since the barriers to the leads are 
assumed to be large, the single-particle wavefunctions of 
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FIG. 1. Conductance vs. chemical potential fj, through a 
linear array of 10 GaAs quantum dots with one spin-1/2 or- 
bital per dot. e 2 /C g = lmeV, d = 0, t — O.lmeV, and 
T — 35mK. Splitting of the two Coulomb blockade peaks 
into minibands is driven by t. The suppression of the 5th 
peak in (b) is the result of a density-dependent SPT. 



the array are like those of a particle in a one-dimensional 
box. The lowest eigenstate has a maximum in the cen- 
ter of the array and a long wavelength, hence a small 
amplitude on the end dots, leading to a suppression of 
the 1st conductance peak. Higher energy single-particle 
states have shorter wavelengths, and hence larger ampli- 
tudes on the end dots, leading to conductance peaks of 
increasing height. The suppression of the conductance 
peaks at the top of the 1st miniband can be understood 
by an analogous argument in terms of many-body eigen- 
states; the 10th electron which enters the array can be 
thought of as filling a single hole in a Mott insulator, etc. 

In Fig. 0(b), the spin-degeneracy of the quantum dot 
orbitals is lifted by the Zeeman splitting. There is a 
critical field B c above which the system is spin-polarized; 
for an infinite ID array with Ci = and t -C £7 fll3|| 
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(27m — sin 27rn), 
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where n < 1 is the filling factor of the lower Hubbard 
band. (Recall that we are here considering only the sin- 
gle spin-1/2 orbital nearest Ep in each quantum dot — the 
magnetic field required to spin-polarize an entire quan- 
tum dot is much larger |H.) Because B c is a function of 
n, one can pass through this spin-polarization transition 
(SPT) by varying n at fixed B. In Fig. |l|(b), this transi- 
tion occurs between the 4th and 5th electrons added to 
the array, consistent with the prediction of Eq. (||). The 
effect of this transition on the conductance spectrum is 
dramatic: The first 4 electrons which enter the array have 
spin aligned with B (up) , but the 5th electron enters with 
the opposite spin, and goes predominantly into the lowest 
single-particle eigenstate for down-spin electrons, which 
couples only weakly to the leads, leading to a suppres- 
sion of the 5th resonant conductance peak by over an 
order of magnitude. It should be emphasized that the 
heights of the conductance peaks change dis continuously 
as a function of B each time there is a spin-flip. 

Splitting of the Coulomb blockade peaks due to inter- 
dot coupling and suppression of the conductance peaks 
at the miniband edges have recently been observed ex- 
perimentally by Waugh et al. H. However, it has been 
pointed out Q that both effects can also be accounted 
for by a model jj| of capacitively coupled dots with com- 
pletely incoherent interdot tunneling. It is therefore of 
interest to consider the effects of interdot capacitive cou- 
pling in the regime of coherent interdot transport. A 
nonzero interdot capacitance Ci introduces long-range 



(1) and decreases 
2 shows the spin 



electron-electron interactions in Eq. 
the intradot charging energy U. Fig. 
susceptibility Xs f° r Ci/C g = 1/2 in linear arrays with 
8 electrons on 12 dots and 10 electrons on 10 dots. The 
n-dependence of B c in Fig. || is qualitatively similar to 
that in a system with intradot interactions only, but the 
values of B c are roughly twice those of a system with 
Ci — 0. Note the rapid growth of \s as B — > B c . In 
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FIG. 2. Spin susceptibility X s = h~ 1 AS/AB at T = vs. 
magnetic field _B for linear arrays of GaAs quantum dots with 
e 2 /C g = lmeV, C t /C g = 0.5, and t = .05meV. Squares: 10 
electrons on 10 dots (B c w 1.5T); triangles: 8 electrons on 12 
dots (B c » IT). 

an infinite array, x s is expected to diverge as £> — > B c 
because the system undergoes a second order quantum 
phase transition Jl3|]. The SPT predicted to occur in 
an array of coupled quantum dots is in contrast to that 
observed in a single quantum dot fL4| , where the critical 
point occurs for minimum total spin. 

Disorder introduces a length scale which cuts off the 
critical behavior as B — > B c . However, as shown in Fig. ^, 
where disorder St ~ t has been included in the hopping 
matrix elements, the SPT has a clear signature in the 
magnetotransport even in a strongly disordered system. 
In Fig. the peak splitting due to capacitive coupling is 
roughly ten times that due to interdot tunneling, so that 
the peak positions are within ~ 10% of those predicted 
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FIG. 3. Conductance vs. chemical potential fi through a 
linear array of 6 GaAs quantum dots with one spin- 1/2 or- 
bital per dot. e 2 /C g = lmeV, Ci/C g = 0.5, i = .05meV, 
T = 120mK. Disorder 5t/i ~ 1 (tit = iix) is present in 
the hopping matrix elements. The splitting of the Coulomb 
blockade peaks into multiplets is dominated by d; however, 
the effect of B is similar to that in Fig. 1. 



FIG. 4. Conductance vs. chemical potential through a 
linear array of 6 GaAs quantum dots with one spin- 1/2 or- 
bital per dot. e 2 /C g = lmeV, d/C g = 0.5, t = .05meV, 
T — 120mK. Spin-dependent disorder St/i ~ 1 (tjf 7^ ti±) is 
included in the hopping matrix elements. Solid curve: B — 0; 
dotted curve: B — 1.3T. At 1.3T, the 2nd conductance peak 
is enhanced by a factor of 1600. 

by a classical charging model ||]. However, the dra- 
matic dependence of peak heights on magnetic field — 
the 4th conductance peak in Fig. ||(b) is suppressed by 
a factor of 32 compared to its B = value due to the 
density-dependent SPT described above — can not be ac- 
counted for in a model which neglects coherent inter- 
dot tunneling. This effect should be observable provided 
g^sB c > max(fcsT, h/r^), where Tj is the inelastic scat- 
tering time. We believe that this striking magnetotrans- 
port effect is the clearest possible signature of a coherent 
molecular wavefunction in an array of quantum dots. 

Fig. ^ shows the conductance spectrum for an array of 
6 quantum dots with the same parameters as in Fig. |^, 
but with spin-dependent disorder in the hopping matrix 
elements, as could be introduced by magnetic impurities. 
Several conductance peaks at B = (solid curve) are 
strongly suppressed due to a many-body enhancement 
of localization. This effect arises because repulsive on- 
site interactions enhance spin-density wave correlations, 
which are pinned by the spin-dependent disorder 
At B = 1.3T (dotted curve) the system is above B c and 
is spin-polarized, circumventing this effect. The second 
conductance peak is enhanced by a factor of 1600 at 1.3T 
compared to its size at B = (not visible on this scale). 
This giant magnetoconductance effect is a many-body ef- 
fect intrinsic to the regime of coherent interdot transport. 

Another interesting phenomenon stemming from the 
competition between coherent interdot tunneling and 
charging effects is the Mott-Hubbard metal-insulator 
transition (MH-MIT), which occurs when collective 
Coulomb blockade (CCB) || is destroyed due to strong 
interdot coupling. For GaAs quantum dots larger than 
about lOOnm in diameter, we find that this transition is 
caused by the divergence of the effective interdot 
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FIG. 5. Conductance vs. chemical potential through a 
linear array of 5 GaAs quantum dots with 5 spin- 1/2 or- 
bitals per dot. e 2 /C 3 = lmeV, A = 0.2 meV, T = .29K, 
Ci n) /Cg = 2 n_1 , and t„ = 0.05 meV(1.05) n (n = 0,...,4). 
The energy gap between Hubbard minibands is not resolved 
for jj, > 9meV (breakdown of CCB). Note the quenching of 
magnetoconductance effects in the ballistic regime. 

capacitance, similar to the breakdown of Coulomb block- 
ade in a single quantum dot |lC( ] . Within the framework 
of the scaling theory of the MH-MIT jjj), one expects 
a crossover from CCB to ballistic transport in a finite 
array of quantum dots when the correlation length £ in 
the CCB phase significantly exceeds the linear dimen- 
sion L of the array. Fig. [| shows the conductance spec- 
trum for 5 quantum dots with 5 spin-1/2 orbitals per 
dot. The divergence of the effective interdot capacitance 
as the interdot barriers become transparent is simulated 
by setting c\ n) jC g = 2"" 1 , n = 0, . . . , 4. In Fig. g, mini- 
bands arising from each orbital are split symmetrically 
into multiplets of 5 peaks by CCB, with the center to cen- 
ter spacing between multiplets equal to e 2 /C g , while the 
energy gap between minibands corresponds to the band 
gap ~ A enhanced by charging effects. The CCB energy 
gap is evident in the first 3 minibands, but is not resolv- 
able for the higher orbitals (Ci/C g > 4), although there 
is still a slight suppression of the conductance peaks near 
the center of the 4th miniband. Comparison of the com- 
pressibility of the system to a universal scaling function 
for the MH-MIT calculated by the method of Rcf. Q in- 
dicates £/L ~ 10 3 for Ci/C g = 8, so that the transport in 
the 5th miniband is effectively ballistic. The peak spac- 
ing within a miniband saturates at e 2 / LC g (plus quan- 
tum corrections ~ t/L) in the ballistic phase because the 
array behaves like one large capacitor, as observed ex- 
perimentally in Ref. Fig. ||(b) shows the effects of a 
magnetic field on the conductance spectrum: a sequence 
of SPTs is evident in the different minibands, with B c 
an increasing function of Ci/C g , leading to quenching of 
magnetoconductance effects in the ballistic regime. 
A finite-size scaling analysis of the compressibility in- 



dicates that the MH-MIT probably occurs at Ci / Cg = oo 
in an infinite array of quantum dots, when the interdot 
barriers become transparent to one transmission mode 

In conclusion, we predict that low-temperature magne- 
totransport experiments on judiciously fabricated quan- 
tum dot arrays would lead to the observation of a vari- 
ety of phenomena resulting from the interplay between 
coherent interdot tunneling and charging effects: SPT, 
MH-MIT, and giant magnetoconductance due to a many- 
body enhancement of localization by magnetic scattering. 
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